Three-dimensional stability of a rectangular plate under uniaxial tension

The paper studies the three-dimensional stability of an isotropic, linear elastic, rectangular plate under a uniform tensile load applied to its sides. The concept of free strains is used to reduce the three-dimensional problem to a two-dimensional one. It is solved using the three-dimensional linearized theory of stability. An approximate solution of the buckling problem is obtained by the finite-difference method. Numerical results are presented __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 7, pp. 116–123, July 2006.     

Interaction of a shock wave with a loose dusty bulk layer

The interaction of a planar shock wave with a loose dusty bulk layer has been investigated both experimentally and numerically. Experiments were conducted in a shock tube. The incident shock wave velocity and particle diameters were measured with the use of pressure transducers and a Malvern particle sizer, respectively. The flow fields, induced by shock waves, of both gas and granular phase were visualized by means of shadowgraphs and pulsed X-ray radiography with trace particles added. In addition, a two-phase model for granular flow presented by Gidaspow is introduced and is extended to describe such a complex phenomenon. Based on the kinetic theory, such a two-phase model has the advantage of being able to clarify many physical concepts, like particulate viscosity, granular conductivity and solid pressure, and deduce the correlative constitutive equations of the solid phase. The AUSM scheme was employed for the numerical calculation. The flow field behind the shock wave was displayed numerically and agrees well with our corresponding experimental results.      

The space of lines in cyclic covers of projective space

Take positive integers m, n and d. Let Y be an m-fold cyclic cover of ${\mathbb{P}}^n$ ramified over a general hypersurface $X\subset P^n$ of degree md. In this paper we study the space $F(Y)$ of lines in Y and show that it is smooth of dimension $2(n-1)-d(m-1)$ if $md>2n-3$ and $2(n-1)-d(m-1)\ge 0$. When $2(n-1)=d(m-1)$, our result gives a formula on the number of m-contact order lines of X (see Definition 1.2).     

Asymptotic error distribution for the Euler scheme with locally Lipschitz coefficients

In traditional works on numerical schemes for solving stochastic differential equations (SDEs), the globally Lipschitz assumption is often assumed to ensure different types of convergence. In practice, this is often too strong a condition. Brownian motion driven SDEs used in applications sometimes have coefficients which are only Lipschitz on compact sets, but the paths of the SDE solutions can be arbitrarily large. In this paper, we prove convergence in probability and a weak convergence result under a less restrictive assumption, that is, locally Lipschitz and with no finite time explosion. We prove if a numerical scheme converges in probability uniformly on any compact time set (UCP) with a certain rate under a global Lipschitz condition, then the UCP with the same rate holds when a globally Lipschitz condition is replaced with a locally Lipschitz plus no finite explosion condition. For the Euler scheme, weak convergence of the error process is also established. The main contribution of this paper is the proof of $\sqrt{n}$ weak convergence of the normalized error process and the limit process is also provided. We further study the boundedness of the second moments of the weak limit process and its running supremum under both global Lipschitz and locally Lipschitz conditions.     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.     

A Runge-Kutta–based WENO reconstruction on moving mesh for compressible fluids

In this paper, we construct a high-order moving mesh method based on total variation diminishing Runge-Kutta and weighted essential nonoscillatory reconstruction for compressible fluid system. Beginning with the integral form of fluid system, we get the semidiscrete system with an arbitrary mesh velocity. We use weighted essential nonoscillatory reconstruction to get the space accuracy on moving meshes, and the time accuracy is obtained by modified Runge-Kutta method; the mesh velocity is determined by moving mesh method. One- and two-dimensional numerical examples are presented to demonstrate the efficient and accurate performance of the scheme.     

A FULLY DISCRETE IMPLICIT-EXPLICIT FINITE ELEMENT METHOD FOR SOLVING THE FITZHUGH-NAGUMO MODEL

This work develops a fully discrete implicit-explicit finite element scheme for a parabolic-ordinary system with a nonlinear reaction term which is known as the FitzHugh-Nagumo model from physiology. The first-order backward Euler discretization for the time derivative, and an implicit-explicit discretization for the nonlinear reaction term are employed for the model, with a simple linearization technique used to make the process of solving equations more efficient. The stability and convergence of the fully discrete implicit-explicit finite element method are proved, which shows that the FitzHugh-Nagumo model is accurately solved and the trajectory of potential transmission is obtained. The numerical results are also reported to verify the convergence results and the stability of the proposed method.     

Stokes equations under nonlinear slip boundary conditions coupled with the heat equation: A priori error analysis

In this work, we consider the heat equation coupled with Stokes equations under threshold type boundary condition. The conditions for existence and uniqueness of the weak solution are made clear. Next we formulate the finite element problem, recall the conditions of its solvability, and study its convergence by making use of Babuska–Brezzi's conditions for mixed problems. Third we formulate an Uzawa's type iterative algorithm that separates the fluid from heat conduction, study its feasibility, and convergence. Finally the theoretical findings are validated by numerical simulations.